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The sky distribution of positronium annihilation continuum 
emission measured with SPI/INTEGRAL* 
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Abstract. We present a measurement of the sky distribution of positronium (Ps) annihilation continuum emission 
obtained with the SPI spectrometer on board ESA's INTEGRAL observatory. The only sky region from which 
significant Ps continuum emission is detected is the Galactic bulge. The Ps continuum emission is circularly 
symmetric about the Galactic centre, with an extension of about 8° FWHM. Within measurement uncertainties, 
the sky distribution of the Ps continuum emission is consistent with that found by us for the 511 keV electron- 
positron annihilation line using SPI. Assuming that 511 keV line and Ps continuum emission follow the same 
spatial distribution, we derive a Ps fraction of 0.92 ± 0.09. These results strengthen our conclusions regarding 
the origin of positrons in our Galaxy based on observations of the 511 keV line. In particular, they suggest that 
the main source of Galactic positrons is associated with an old stellar population, such as Type la supernovae, 
classical novae, or low-mass X-ray binaries. Light dark matter is a possible alternative source of positrons. 
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>• ! 1. Introduction 



The annihilation of positrons with electrons gives rise 
to two characteristic emissions at gamma-ray ener- 
gies: the hallmark line at 511 keV, and the unique 
three-photon positronium (Ps) c ontinuum emission (cf. 
iGuessoum. Jean, fc Gillardll2005|) . Direct annihilation of 
positrons with electrons, and their annihilation via the for- 
mation of para-Ps (with the spins of electron and positron 
being anti-parallel), result in the emission of two 511 keV 
photons. Annihilation via the formation of ortho-Ps (with 
the spins of electron and positron being parallel) produces 
three photons and gives rise to the Ps continuum emis- 
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sion, which is roughly saw-tooth shaped with a peak at 
the maximum energy of 511 keV l|Ore & Powelllll94^) . 

Cosmic positron annihilation radiation was first de- 
tected from the Galactic centre (GC) direction in bal- 
loon observations during the 1970s and has been the 
focus of intense scrutiny by a large number of balloon 
and satellit e borne expe r iments ever since (see e.g. the 
reviews bv lTueller| Il992t iHarrisl HqQtI . Despite tremen- 
dous observational and theoretical efforts, the origin of 
the positrons is still poorly understood. A large variety 
of positron sources and produc tion mechanisms have been 
proposed over the years (e.g. IChan fc Lingenferterlll993t 
iDermer fc Murphvl 12001^1 . Among the more promising 
source candidates are radioactive nucleosynthesis prod- 
ucts from supernovae. More recently , hypernovae/GRBs 
(Cassc et al. 2004; Par izot et al . 2005) and light dark mat- 
ter fe.g. lBoehm et al.ll2004j) have been proposed as possi- 
ble candidates. Compact objects comprise another poten- 
tial candidate source of positrons. 

Investigations of the sky distribution of the annihi- 
lation radiation promise to provide clues to the identi- 
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fication of the source(s) of positrons in our Galaxy, de- 
spite the fact that positrons may travel from their birth 
places before annihilating. First maps of the annihila- 
tion radiation, limited to the inner regions of our Galaxy, 
were obtained using the OSSE instrum ent on board 
the C ompton Gamma-Ray Observatory l|johnson et alJ 
Il993l) in the 511 keV line and in Ps continuum emis- 
sion ( e .g. IPurcel et alll997t IChen et alll997t iMilne et alJ 



1200(1 l2nnialbl. I20r)2ft . Fu r thermore, t he OSSE instru- 
ment allowed iKinzer et"al1 ( 19991 l200ll) to study the one- 
dimensional distribution in longitude and in latitude of 
diffuse emission, including annihilation radiation, from 
the inner Galaxy. With the commissioning of the imag- 
ing spectrometer SPI on board ESA's INTEGRAL obser- 
vatory, high spectral resolution mappi ng with improved 
angular resolution has become feasible pe an et al.|| 2003at 
IWeidenspointner et alJl2004l IKnodlseder et al l2005|l -The 
511 keV line emission is found to be dominated by the 
Galactic bulge and/or halo; emission from the Galactic 
disk is much fainter, implying that positron annihila- 
tion is concentrated in the central regions of our Galaxy 
iKnodlseder et alJl2005lh 

A first SPI measurement of the flux in Ps contin- 
uum emission, and of the Ps fraction fp s (the frac- 
tion of positrons that annihilate through Ps forma- 
tion: iBrown & Leventhall Il987h . has been presented by 
[Ch urazov et alJ (2005). The value of fp s , as well as the 
detailed shape of the 511 keV annihilation line, depend 
on the physical properties of the annihilation media; 
therefore detailed spectroscopy of the positron annihila- 
tion can provide uniq ue information on the annih i lation 
media and proc esses i Guessoum. Jean, fc Gillardl 12005: 
I Jean et al-lEoO^ . 

In this publication, we present results concerning the 
Galactic distribution of Ps continuum emission using ob- 
servations of most of the celestial sphere with the spec- 
trometer SPI/INTEGRAL. The analysis of the SPI obser- 
vations is presented in Sec. [3 O ur mapping, model fitting, 
and spectral results are given in Sec. [3] A summary and 
our conclusions can be found in Sec. 



2. Instrument description and data analysis 

The SPI imaging spectrometer consists of an array of 19 
actively cooled high resolution Ge detectors, which cover 
an energy range of 20-8000 keV with an energy resolution 
of about 2.1 keV FHWM at 511 keV. SPI employs an 
active anti-coincidence shield made of bismuth germanate 
(BGO), which also acts as a collimator. In addition to 
its spectroscopic capabilities, SPI can image the sky with 
moderate spatial resolution of about 3° FWHM using a 
tungsten coded aperture mask. The fully coded field-of- 
view of the instrument is about 16° . A detailed descrip tion 
of the instrument was given bv lVedrenne et alJ l)2003|) . 

The analysis presented here is based on the same 
data set that we used to study the sky distribution of 
the 511 keV electron-positron a nnihilation line (for de- 
tails see IKnodlseder et a.1.1 l2005h . The data comprise all 



observations included in the December 10, 2004 public 
INTEGRAL data release (i.e. three-day orbital revolu- 
tions 16-76, 79-80, 89-122) supplemented by instrument 
team observations of the Vela region during revolutions 
81-88. The observations were taken during the epoch 
December 9, 2002 through September 12, 2003. In total, 
the data set consists of 6821 pointings 1 with a combined 
live time of 1 3.5 x 10 6 s. The resultin g exposure to the sky 
(see Fig. 1 in IKnodlseder et aI1l2004 is relatively uniform 
for Galactic longitudes \l\ < 50° and Galactic latitudes 
\b\ < 15°, with a maximum at the GC region. Unexposed 
regions arc limited to the south Galactic pole region and 
a few areas at intermediate Galactic latitudes. 

The analysis of the Ps continuum emission, and of 
the underlying Galactic power-law continuum (hereafter: 
Galactic continuum) emission 2 , is rendered challenging by 
a signal-to-background ratio which is only of the order of 
a few tenths of one per cent (compared to a few per cent 
for the 511 keV annihilation line). Background modelling 
therefore is crucial. Our approach consists of fitting time 
series of detector count rates by a linear combination of 
so-called background templates (time series of background 
components, explained below). The detector signal arising 
from a celestial source is also variable in time, even if the 
source is intrinsically stable, because SPI's exposure to it 
varies in time as the instrument performs a series of obser- 
vations. Source components can be distinguished from the 
dominant background components by taking advantage of 
their differing time series. 

Many processes contribute t o the ins t rumen- 
tal background at all e nergies l|Jean et alJ l2003bt 
IWeidenspointner et alJl2003|h The statistical significance 
of the available data is still a limiting factor in the 
analysis of sky signa ls. We found that the b ackground 
components used by IKnodlseder et al.l l|2005h to model 
the instrumental 511 keV line provide a viable compro- 
mise between accuracy and simplicity. The model consists 
of three components: a constant, the rate of saturating 
(> 8 MeV) events in the Ge detectors, and a template 
representing the exponential build-up of long-lived iso- 
topes (a lifetime of 352 days, corresponding to 65 Zn, was 
chosen 3 ). The role of the first two templates is to account 
for prompt and short-lived background components, 



1 Each INTEGRAL observation is split into numerous so- 
called pointings, often arranged on a square grid of 25 posi- 
tions on the sky centred on the main target of the observation, 
which each pointing lasting about 2000 s iCourvoisier et aT] 
20031). This so-called dithering strategy is necessary to mini- 
mize imaging artefacts for SPI by providing a sufficient number 
of measurements (sky signal and instrumental background) for 
each of the 19 detectors to constrain sufficiently the fluxes from 
each pixel in the field-of-view. 

2 This emission is partly due to cosmic-ray interactions in 
the interstellar medium, partly due to point sources. 

3 Since our data extend only over about 1 year, we cannot 
yet distinguish different half-lives longer than about 6 months. 
The exact value of the lifetime is therefore not crucial, but it is 
important to allow for the significant long-term build-up effect. 
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the third template accounts for long-lived background 
components. This three-component background model 
was applied either separately to three energy intervals, 
namely 410-430, 447-465, and 490-500 keV (hereafter: 
Ps analysis intervals), or to the summation of the events 
in these three intervals into a single data set (hereafter: 
summed Ps analysis intervals). The intervals were chosen 
to be free of strong background features (which require 
specific modelling using other, intermediate time scales). 
We do not include the fourth background component 
used in the 511 keV line analysis, which represented the 
continuum underlying the instrumental 511 keV line. 
Thus we are measuring the total emission, Ps continuum 
plus Galactic continuum. 

As in studies of the sky distribution of the 511 keV 
line emission, we found that by fitting the normaliza- 
tions of each of the three background model compo- 
nents for each detector remaining systematic uncertain- 
ties introduce a bias. However, if we fit the second 
background component not only for each SPI detector, 
but additional ly for each orbital revo lution (as in model 
0RBIT-DETE in IKnodlseder et al.ll2005l) we reduce system- 
atic trends to an acceptable level. This method is sim- 
ilar to those appli ed by us in earlier analyses of the 
511 keV line (e.g. IJean et all l2003al IKnodlseder et alJ 
l2003HWeidenspointner et al.ll2004h ~ 

To obtain model independent maps of extended emis- 
sion from our Galaxy we employed an implementation of 
th e Richardson-Lucy algor ithm that is described in detail 
in IKnodlseder et all {2005) . Briefly, in this approach the 
intensity attributed to each sky pixel is calculated itera- 
tively using the relation 

where k is the iteration number, fj is the intensity from 
sky pixel j in iteration k, Rij is the instrumental re- 
sponse matrix (linking data space bin i to sky pixel j), 
rii is the number of counts measured in data space bin i, 
ef = J2jLi RijJj + bi is the predicted number of counts 
in data space bin i after iteration k (with bi being the 
predicted number of instrumental background counts for 
data space bin i), N and M are the dimensions of the data 
and image spaces, respectively, and X k is an acceleration 
factor. To avoid noise artefacts in weakly exposed regions 
of the sky, the image increments are weighted with a quan- 
tity uij = (X^iLi Rij) 2 that is related to the sensitivity of 
the instrument. 

When fitting models of intensity distributions on the 
sky to the data, we employed the maximum likelihood 
multi- component model fitting algorithm described in de- 
tail in IKnodlseder et aD (|^005). Briefly, assuming Poisson 
noise for the measured number rii of events in each of 



the N data space bins, the algorithm maximizes the log 
likelihood function 

N 

In L = rii In £i — e.i — In rii}. (2) 

i=l 

where = J2k a k s i + bi(f3) is the predicted number of 
(source plus background) counts in data space bin i, sf = 
J2j=i fjRij m * ne s ^y intensity model k folded into data 
space bin i (with /j° being the intensity in sky pixel j due 
to sky model k and Rij being the instrumental response 
matrix), 6;(/3) is the background model, and ak and (3 are 
scaling factors for the sky intensities and the background 
model components, respectively, that are adjusted by the 
fit. 

3. Results 

3.1. Mapping 

We employed an implementation of the Richardson-Lucy 
algorithm, outlined in Sec. [3J to obtain a model indepen- 
dent map of the extended emission in the summed Ps 
analysis intervals (i.e. the events were first summed and 
then analyzed, hence the total emission in the three en- 
ergy intervals is mapped without constraints on the rel- 
ative fluxes in each of them; spectroscopy of the emis- 
sion is presented in Sec. l3.3[l . Since we are only interested 
in extended emission, the two point sources Crab and 
Cygnus X-l (which are clearly seen due to their strong 
continuum emission at these energies) were subtracted in 
the image reconstruction process. In addition, to reduce 
noise artefacts, the iterative corrections were smoothed 
during image reconstruction with a 5° x 5° boxcar aver- 
age. The resulting sky map after iteration 8 is presented 
in Fig. n We selected this iteration because at this point 
of the iterative mapping procedure the fit quality corre- 
sponds approximately to that which we obtained when 
fitting sky models to the data (see Sec. I3.2fl . This choice 
ensures that the algorithm has not yet overfitted the data, 
which could result in substantial image noise and image 
artefacts. We note, however, that analyses of simulated 
observations show that at this point in the iterations and 
with the limited amount of data available, features such 
as the faint extended emission expected from the Galactic 
plane would not necessarily be apparent, and indeed the 
Galactic plane is not apparent in Fig.Q]although we detect 
it at a low significance level in model fits (see Sec. l3.2JI . The 
simulations also show that with increasing iterations (in 
our analysis beyond iteration 8) faint extended emission 
begins to appear, however, so do artefacts due to overfit- 
ting, and a distinction in the image between actual signal 
and noise becomes increasingly difficult. 

The only prominent signal seen is that from the 
Galactic bulge region. The surface brightness of any emis- 
sion from any other sky regions is much fainter. The emis- 
sion appears to be symmetric about the GC, and its cen- 
troid coincides well with the GC. In short, the character- 
istics of the extended emission in the three Ps analysis 
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Fig. 1. A Richardson-Lucy sky map of extended emission in the summed Ps analysis intervals (the combination of 
the intervals 410-430, 447-465, and 490-500 keV). The contour levels indicate intensity levels of 10~ 2 , 10~ 3 , and 10~ 4 
ph cm -2 s _1 sr _1 . Details are given in the text. 



intervals agree well with those already identified by us 
for the 511 keV line radiation. This agreement suggests 
that the extended emission from the GC region in the 
Ps analysis intervals is dominated by Ps continuum emis- 
sion, with the Galactic continuum emission being fainter 
(as will be demonstrated in Sec. 13. 3[) . We note that these 
characteristics are consistent with OSSE maps obtained 
for th e Ps continuum emission fe.g. IMilne et al.ll2001bl 
in which the surface brightness is highest in the 

bulge. 

We also used the SPIROS algorithm 
l|Skinner fc Connelll Effll to search for point sources in 
these three energy bands. No significant sources in addi- 
tion to the Crab and Cygnus X-l were found outside the 
GC region. When imaging the GC region with SPIROS, 
the observed emission formally is equally well explained 
by six point sources. However, different sets of sources are 
found in the different Ps analysis intervals. These sources 
do not coincide with any known high-energy sources, 
including the four h ighest energy sources reported by 
iBouchet et, alJ (|2005j) . The spectrum of these four sources 
drops steeply above about 300 keV, and since we are 
analyzing rather narrow energy intervals above 400 keV 
the fact that we do not yet detect them is not surprising. 
We therefore conclude that the point sources found by us 
using SPIROS are all spurious, resulting from SPIROS' 
attempt to account for intrinsically diffuse emission with 
a set of point sources. 



3.2. Model fitting 

A more quantitative approach for studying the Galactic 
distribution of the observed extended emission is model 
fitting, which we performed using a maximum likelihood 



multi -component fitting algorithm l|Knodlseder et all 
l2005h outlined in Sec. El 

We first modelled the emission in the three summed 
Ps analysis intervals 4 by an ellipsoidal distribution with a 
Gaussian radial profile and determined the best-fit cen- 
troid location (lo,b ) an d extent in Galactic longitude 
and latitude (FWHM;, FWHM b ). We then combined this 
Galactic bulge model with one of t wo models for emis- 
sion from the G alactic disk: both HI i Dickev fc Lockmanl 
Il990l) and CO l|Dame et al.lll987f) distributions are trac- 
ers of Galactic mat ter and are believed to correlate with 
diffuse emission fcf. lHarris et al.lll990UKinzer et alJ ll999: 
Str ong et alJ I2004j) . The results of these fits are sum- 
marized in Table ^ In each of these fits, the Crab and 
Cygnus X-l were included as steady point sources whose 
intensities were fitted. W hen including the fou r highest- 
energy sources reported bv lBouchet et al.l l|2005|) the qual- 
ity of the fits is only slightly improved and the fit results 
do not change significantly; therefore these point sources 
were excluded from the final analysis. 

As can be seen from Tabled the centroid of the bulge 
emission is the same within errors for all three models. 
There is marginal evidence for a slight offset of the cen- 
troid from the GC, but it is of a magnitude that could 
easily result from the combined effects of statistical and 
systematic biases in the background model (indeed, there 
is a similarly marginal, but o pposite, offset of the cen troid 
in the 511 keV line emission; iKnodlseder et al.l2005^ . The 
extent of the bulge emission, and its flux, do depend on the 
sky model. If the extended emission is modelled by a bulge 
component only, then there is marginal evidence for the 
bulge emission to be more extended in longitude than in 



4 Results for the individual energy intervals are consistent 
within statistical uncertainties. 
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Table 1. A summary of results from model fitting in the 
summed Ps analysis intervals (the combination of the in- 
tervals 410-430, 447-465, and 490-500 keV). In the table, 
A is the maximum likelihood ratio of the fit, and 
fco are the fluxes in the bulge and the two different disk 
components, respectively. The remaining symbols are de- 
fined in the text. The error bars reflect the statistical Ict 
confidence interval for one degree of freedom. 



Quantity 



Value 



Bulge only 



A 

io n 

bo [°] 

FWHMi 

FWHMb 



494.5 
1.3 ±0.8 
-0.8 ±0.6 



10.0 



+4.1 



f b [ph/cm7s] 1.20l°;i* x 10 



6.6ii:* 

+0.22 
0.20 

+0.14 
0.11 



0.66Z 



Bulge ± HI 



A 

lo [°] 
bo [°] 
FWHM, 
FWHM 6 
c 



502.2 
1.3 ±0.8 
-0.9 ±0.6 

s.3±!:S 

6.8l 15 
0.7lj 



1.2 
+0.26 
0.20 

16±S"12 x 10" 



f b [ph/cm /s] ^"_ . 10 
fm [ph/cm 2 /s] 4.19t 1 1 j 2 x 10 



Bulge + CO 



A 

i n 
^ n 

FWHMi 
FWHM 6 
e 

f b [ph/cm 2 /s 
fco [ph/ 



510.5 
1.3 ± 1.1 

-1.2 ±1.0 

92 +3.i 

7.. 



cm 2 /s] 



^-2.1 
j+3.2 

0.85+°;^ 
0.86^'^ x 10~ :i 
1.92tgig x IO" 3 



latitude (the ellipticity e = FWHM b /FWHM ; deviates by 
about 1.5cr from unity). However, inclusion of a Galactic 
disk component improves the fits, with the significances of 
the HI distribution and of the CO distribution being about 
2.8(7 and 4. Oct, respectively, favouring the latter. Another 
reason to adopt the CO distribution as the better disk 
model of the two is the fact that the resulting total sky flux 
of about (2.8 ± 0.5) x 10~ 3 ph cm~ 2 s _1 agrees well with 
the value of about 2.5 x 10 -3 ph cnr 2 s^ 1 determined with 
SMM 5 in the Ps analysis intervals, whereas the total bulge 
and HI disk model flux of (5.4 ± 1.5) x 10 -3 ph cm~ 2 s _1 
is only marg i nally consistent with the SMM spectrum of 
lHarris et aill|l990lh 

Inclusion of a Galactic disk component in the fits also 
renders evidence for ellipticity of the bulge component 

5 The Gamma Ray Spectrometer on board the Solar 
Maximum Mission jForrest et alJll986Tl . 



insignificant. The bulge shape is consistent with circular 
symmetry, with a FWHM of abo ut 8°, in agreement wit h 
our results for the 511 keV line (K nodlseder et alJl2005h . 
As is the case for the annihilation line, the extent of the 
Ps continuu m bulge emission i s slightly larger than that 
derived bv iKinzer et alJ l|200lj) from OSSE observations. 
However, the difference is not very significant, and it is 
possible that there is b ias in the OSSE an alysis favouring 
a smaller bulge extent l|Kinzer et alJl200l|) . 

The fluxes that are attributed to the disk components 
exceed the bulge flux by factors of 2-4 (see Table P). 
However, since the disk flux is distributed over a much 
larger sky region, the corresponding surface brightness is 
much lower. The model fits therefore confirm the mapping 
result: the intensity of extented Galactic emission in the 
three Ps analysis intervals is brightest in a region, several 
degrees across, around the GC. 

3.3. Spectral analysis 

To assess the relative contributions of Ps continuum and 
Galactic continuum emission to the emission mapped in 
Fig.^ we performed a spectral analysis in six energy inter- 
vals: the three bands used for the Ps continuum analysis 
supplemented by the energy intervals 500-507, 507-515, 
and 515-565 keV. More det ailed spectroscopy including 
the line shape is presented in I Jean et alJ l|2006|h We want 
to emphasize that present uncertainties in the sky distri- 
bution of the emission components, in particular for the 
Galactic continuum emission, imply systematic uncertain- 
ties for this kind of spectral analysis. Similarly, when com- 
paring different results, one has to keep in mind that there 
will be a systematic difference between fluxes obtained un- 
der different assumptions for the spatial distribution of the 
emission. 



3.3.1. Simplistic SPI spectroscopy 

A first, and rather simplistic, spectral analysis of the 
SPI result for the diffuse emission from the GC region 
(I^Ij \b\ < 16°) was performed by fitting three spectral com- 
ponents: a line at 511 keV, the Ps continuum, and a power 
law 6 representing the Galactic continuum. Before fitting, 
these models were convolved with the SPI instrumental 
resolution. Following our best sky model ("Bulge±CO", 
see Table ^1 , and in agreement with our results on the 
sky distribution of the 511 keV line emission, the SPI 
fluxes used in this fit were obtained assuming that the 
bulge emission is circularly symmetric about the GC with 
a Gaussian radial profile and a FWHM of 8°, and assum- 
ing that the disk emission is traced by the line-of-sight 
integrated CO distribution (hereafter: model G8CO). 

The limited data available for this spectral analysis re- 
quired that some parameters of our spectral model be fixed 
or constrained. The ccntroid and width of the astrophys- 



The power law is defined as 
f Y {E) = A ■ (_E/300 keV)- a 
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Fig. 2. A fit of the SPI result for the diffuse emission 
from the GC region (|Z|, |6| < 16°) obtained with a spatial 
model consisting of an 8° FWHM Gaussian bulge and a 
CO disk. In the fit a diagonal response was assumed. The 
spectral components are: 511 keV line (dotted), Ps contin- 
uum (dashes), and power-law continuum (dash-dots). The 
summed models are indicated by the solid line. Details of 
the fitting procedure are given in the text. 

ical positron annihila tion line were fixed at 511 keV and 
2.5 k eV FWHM (c.f. iLoniou et al1l2004t IChurazov et"aD 
2005|, respectively. In addition, t he power-law pho- 
ton index a was fixed at 1.75 (c.f. iKinzer et all Il999t 
iBouchet et al.ll2~005|) . In this first analysis, off-diagonal el- 
ements in the SPI response, which are small, were ne- 
glected. When fitting the spectral models to the data, the 
models were averaged over the energy intervals pertaining 
to each data point. A rigorous spectral analysis, using the 
full response matrix of SPI, will be presented below. The 
advantage of this first analysis is its simplicity, and the 
fact that its results can easily be displayed graphically. 

The best fit result is shown in Fig.0 We find a 511 keV 
line flux of (1.02 ± 0.10) x 10~ 3 ph cm" 2 s - \ and a Ps 
continuum flux of about (5.41 ±1.59) x 10 -3 ph cm~ 2 s^ 1 . 
From these fluxes, and using simple error propagation, we 
obtain a v alue of fp R of 1.04 ± .07 fo r the Ps fraction as 
defined by iBrown fc Leventhall l)l987h . The power law is 
not significantly detected with a normalization of (1.47 ± 
0.93) x 10~ 5 ph cm- 2 s" 1 keV" 1 . The X 2 of the fit is 
about 1.8 for 3 degrees of freedom. Based on the best fit 
parameter values, the Ps continuum emission accounts for 
about 71% of the total emission in the three Ps continuum 
analysis intervals mapped in Fig. ^ 

3.3.2. Rigorous SPI spectroscopy 

To perform a rigorous spectral analysis of the emission 
from the Galactic bulge, we fitted the same spectral model 
- a line at 511 keV, the Ps continuum, and a power 
law - to a grid of points centred on the GC using the 
INTEGRAL/SPI specific capa bilities of the XSPEC (ve r 
12) analysis package (see e.g. iDorman fc Arnaudll200ll) . 
which allowed us to take into account the full spectral 



and spatial response of the instrument. For a detailed de- 
scription of t h is me thodology the reader is referred to 
IShrader et alJ l|2004h : here we will just briefly describe 
the basic procedure followed. This is the first time that 
XSPEC has been applied to spectroscopy of an extended 
sky source observed with the SPI instrument. As an aside 
note, we wish to warn the interested reader that we found 
the original Ps continuum model in XSPEC, P0SM, to be 
incorrectly implement ed. We developed an d tested a new 
implementation of the lOre fc PoweliTl| 19491) spectral shape 
of Ps continuum emission, which will be included in sub- 
sequent releases of XSPEC. 

The data selected for this portion of our analysis com- 
prise a subset of the total data presented in this paper. 
Observations were selected for inclusion in our spectral fit- 
ting when the SPI telescope axis was aligned with the GC 
to within an angular offset of 16° (the extent of the nomi- 
nal fully-coded SPI field-of-view) . This resulted in a total 
of about 750 spacecraft pointings (Science Windows), to- 
talling ~1.7 Ms of live time, being used in this analysis. 

The full SPI instrument response, including diago- 
nal plus off-diagonal matrix elements, was then com- 
puted, accord i ng to the methodologies described in 
ISturner et al for each SPI detector for each 

selected instrument pointing for each of our grid 
points spatially sampling the bulge region. Specifically, 
we computed the response for a 21-point raster at 
(l,b) = (0 ,0 ),(±4°,0 ),(0 ,±4°),(±8°,0 ),(0°,±8 ), 
(±4°, ±4°), (±8°, ±4°), (±4°, ±8°). 

The data were then simultaneously fitted to the phys- 
ical model described above - 511 keV line, Ps continuum, 
and power law - and the 3-component background model 
described in Sec. [21 The background model in this case 
was parameterized so that small (±10%) variations were 
allowed for the normalization terms of each component in 
each energy interval, using the results of model fits (as de- 
cribed in Sec. 13.2(1 to initialize the background model pa- 
rameters. In practice we found that the background mod- 
elling worked quite well, with the best fit solutions typ- 
ically corresponding to normalization terms within ±1% 
of unity. 

We then made the assumption that the net flux con- 
sists of additive contributions from the two spatial models 
discussed in Sec. l3.3.H i.e. the Gaussian and CO distribu- 
tions of spatial model G8CO. The spectral model was then 
applied to the SPI instrument response function twice at 
each spatial raster, with a normalized, relative, weighting 
factor based on both the Gaussian and the CO distribu- 
tions. This leads to a data space which scales as: (num- 
ber of SPI pointings) x (number of detectors) x (number 
of spectral channels). This number is then multiplied by 
(number of spatial rasters) x (2 spatial distribution mod- 
els) to give the number of individual response matrices ap- 
plied to the spectral model for the x 2 minimization prob- 
lem. This leads to - 750 x 19 x 6 x 21 x 2 - 3.6 x 10 6 folded- 
model calculations per iteration step of the x 2 minimiza- 
tion procedure. Specifically, we used the XSPEC "FLUX" 
command and the best fit parameters of each individual 
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model component to integrate over the covered energy 
range. 

The parameter space was constrained as follows. The 
centroid and width of the positron annihilation line were 
fixed at 511 keV and 2.5 keV FWHM, respectively, as in 
our first analysis (see Sec. 13.3. Q . We fixed the power-law 
photon index a to a value of 1.75, but allowed the am- 
plitude to vary by about a factor of 4 relative to that 
obtained in our first analysis described above. Otherwise, 
the model parameters - specifically the Ps continuum and 
Gaussian line normalization terms - were allowed to vary 
freely in the \ 2 minimization. These two normalization 
terms were varied separately with respect to the two spa- 
tial distributions, but linked from grid point to grid point 
within a given spatial model. This leads to 6 free physi- 
cal model parameters (3 normalizations for each of the 2 
spatial models), in addition to the 18 background model 
parameters (3 parameters in each of the 6 energy intervals) 
for the overall fit. 

We obtained a Ps continuum normalization of (3.11 ± 
0.56) x 1CP 3 ph cm~ 2 s _1 . Combined with the in- 
ferred Gaussian line component normalization of (9.35 ± 
0.54) x 10~ 4 ph cm~ 2 s _1 we obtain a Ps fraction of 
fp s — 0.92 ± 0.09. The normalization of the power-law 
component, rescaling the XSPEC result to the power- 
law function defined in footnote 6, is (3. 79^^25) X 
10~ 6 ph cm~ 2 s _1 keV -1 . Thus over the range of our 
data, the power-law component contributes approximately 
6% (and possibly as much as 14%) of the continuum flux 
based on our model fitting. If we allow the power-law in- 
dex to vary freely, the resulting power-law flux remains 
within the confidence interval above; hence our conclusion 
regarding the flux contribution of the power-law compo- 
nent is robust. The background normalization terms, as 
noted, were within 1% of unity. The \1 value obtained 
was of order unity; specifically, using the full 1.65 Ms 
of the data selected for this analysis, a \ 2 P er degree 
of freedom of 99065.1/86289 w 1.15 was achieved. The 
uncertainties for a given parameter, specficially the line 
and Ps continuum fluxes and the power-law normaliza- 
tion, were derived by varying the parameter within its 
allowed range. At each step, the other free parameters 
are allowed to vary until the fit statistic is minimized, 
determining the ler confidence region for each parame- 
ter (specifically, this is accomplished using the "ERROR" 
procedure of XSPEC vl2). We note that the uncertainty 
in the Ps fraction includes both the variances and the co- 
variances of the 511 keV line and Ps continuum fluxes in 
the variance-covariance matrix of the fit. 

Clearly, given the nature of the data and parameter 
spaces of the problem, it is difficult to present a graphi- 
cal representation of the results. A procedure commonly 
used within XSPEC is to "group" the data for plotting 
purposes, which produces what is essentially an averaged 
representation of the data and model convolved with the 
response. In this application, however, this is not a vi- 
able option as we would be averaging over some ~ 10 6 



spectra/responses/model calculations. We have thus cho- 
sen not to include a plot of the XSPEC results. 

It follows from our spectral analysis of the Galactic 
bulge emission that at least 76% (the best value is 85%), of 
the observed flux in the summed Ps analysis intervals can 
be attributed to Ps continuum emission - which therefore 
dominates the emission from the GC region. The map in 
Fig. n is dominated by emission from the central regions 
of our Galaxy. Hence we can conclude that Fig. ^ provides 
a view of the sky distribution of Ps continuum emission. 



3.3.3. Comparison to SPI results 

We note that our 511 keV line flux is somewhat 
lower than the previous determination of (1.05 ± 
0.06) x 10~ 3 ph cm~ 2 s" 1 for the Galactic bulge by 
iKnodlseder et all l|2005h at the ~ 2a level (the uncer- 
tainties of other ea rlie r SPI measurements, e . g. tho se by 
IJean et all l)2003a|) or IWeidenspointner et al.l (|200^ . are 
too large to reach a definitive conclusion), and that the 
511 keV line, the Ps continuum, and the Galactic contin- 
uum fluxes from our XSPEC analysis are lower than those 
obtained in our more simplistic approach described above. 
We suggest that this may be due to the fact that our 
XSPEC analysis included the off-diagonal response matrix 
elements, while the previous analyses used a diagonal ap- 
proximation. Roughly speaking, the reason for a lower flux 
in our case is that with more matrix elements a smaller 
model flux is required to match the "model counts" with 
the actual detector counts in the \ 2 minimization proce- 
dure. We find for example, that for a Crab-like power-law 
model applied to the SPI full response leads to a ~ 8% 
lower flux than for the same model applied to a diag- 
onalized version of the response. Indeed, for the actual 
Crab data used in determining the SPI calibration, the dis- 
crepancy between flux dete rminations made with SPI ROS 
(diagonal approximation; ISkinner fc Connelll l2003|) and 
XSPEC was initially about 15% and in a similar sense 
(XSPEC flux was lower). Subsequent efforts were made 
to calibrate the S PIROS method for better agreement 
l|Sizun et alJ 20041 . Finally, we note that the effects of the 
off-diagonal terms are even more pronounced for the case 
of a positively sloped spectrum, such as, notably, the ex- 
pected Ps continuum form. 

Our result for the Ps fraction of f Ps = 0.92 ± 0.09 
agrees very well with an ear l ier res ult of fp s — 0.96 ±0.06 
derived bv IChurazov et all l)2005h for the Galactic bulge 
using SPI data. 

We also note that lStrong et al.1 l(2005h concluded, as we 
do, that the emission from the Galactic bulge is dominated 
by Ps continuum emission. Their results for \l\, \b\ < 10° 
imply that about 81% of the emission in the Ps analysis 
intervals is due to Ps continuum emission, in excellent 
agreement to our result. 
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Fig. 3. A fit of the GC s pectrum obtained from OSSE data 
bv IKinzer et alJ l|l999h to the SPI result for the diffuse 
emission from the GC region (|Z|, |6| < 16°) obtained with 
a spatial model consisting of an 8° FWHM Gaussian bulge 
and a CO disk. Details are given in the text. 



3.3.4. Comparison to OSSE results 

The comparison of the astrophysical emission models de- 
termined with other instruments with our SPI results was 
performed by first convolving the spectral models with 
the SPI instrumental resolution, and then by averaging 
the convolved models over the energy ranges of the SPI 
data points. As in our spectroscopic analyses described 
above, the SPI results used in these comparisons were 
obta ined with spatial m odel G8CO. It is worth noting 
that lHarris et al' had assumed in their SMM anal- 

ysis that the Galactic gam ma-ray emission follo ws the CO 
line-of-sight distrib ution of iDame et al.l l|l987|) . and that 
IKinzer et alJ 1 19991) had concluded from their OSSE anal- 
ysis that the longitude distribution of the Galactic con- 
tinuum emission follows that of CO, albeit with a some- 
what broader distribution in latitude. To make compar- 
isons with OSSE results, we consider separately the form 
of the spectrum, the annihilation flux, and the Galactic 
continuum emission. 

A comparison of the Galactic spectrum seen by SPI 
in a box define d by \l\Ab\ < 16° for our spatial model 
G8CO with the lKinzer et all l|l999h OSSE best fit model 
is depicted in Fig.|3 In this comparison of spectral forms, 
the normalization of the OSSE model was obtained by a fit 
to the SPI data. We find a normalization factor of about 
0.6, and a % 2 value of about 3.0 for 5 degrees of freedom. 
We conclude that the shape of the total SPI spectrum 
agrees well with that of the total OSSE spectrum for the 
GC region. 

For comparison with the annihilation radiation anal- 
ysis in the c entra l square radian around the GC by 
IKinzer et all l)200ll) using OSSE data, we repeated the 



simplistic spectral analysis of Sec. I5.3.1 | in \l\,\b\ < 29° 
using spatial model G8CO. IKinzer et alJ l|200l|) obtained 
total annihilation radiation fluxes (i.e. the sum of the 
annihilation line and the Ps continuum fluxes) of about 



(6 — 11) x 10~ 3 ph cm~ 2 s _1 when fitting nova ellipsoid or 
Gaussian bulge models combined with different thin disk 
models, which, of all their models, are the most similar to 
our G8CO model. The SPI total annihilation flux is about 
7.5 x 10 -3 ph cm~ 2 s _1 , in very good agreement with their 
results. Our measurement of the Ps fraction is also in very 
good agreement with their value for the Galactic bulge of 
f Ps = 0.93 ± 0.04. 

For the Galactic continuum emission comparison we 
computed the flux in the 410-500 keV band, that is 
mos t important i n our analysis. C omparison with both 
the IKinzer et ail l|200l|) and the IPurcell et all l|l996ft 
analyses of OSSE observations shows good agreement. 
App roximating the cos mic-ray interaction model used in 



Approximating t he cos mic-ray interaction model used n 
the lKinzer et al.l ( 20011) analysis by a power law with pho 



ton index 1.75, we find that their results correspond to a 
flux of about 3.4 x 10~ 4 ph cm 2 s 1 rad 1 , whereas as- 
suming that our results refer to a 32° range in Galactic 
long itude we find 3.0 + i 3 x 10~ 4 ph cm~ 2 s _1 rad -1 . 
The IPurcell et alJ l|l99fif) analysis implies a flux of about 
3.6 x 10 -4 ph cm -2 s" 1 rad" 1 , however, their power-law 
index of 2.5 is steeper than our result, which we think is 
due to the fact that they used data dominated by energies 
below a few hundred keV. 

3.3.5. Comparison to SMM results 

A comparion of SPI results for the diffuse emission from 
our Galaxy obtained with spatial model G8CO with re- 
sults for the total emissi on from the i nner Galaxy ob- 
tained from SMM data bv lHarris et al ] l|l99ri() (which arc 
not essentially different from the emission from the whole 
Galaxy) is given in Fig.0] As can be seen, the SPI results 
agree well with the best fit SMM model. Except for the 
511 keV line interval o ur fluxes agree with in their statisti- 
cal lcr errors with the lHarris et al ] lll990l) best fit model. 
In the 511 keV line region our flux falls about 2a short 
of the SMM result. Formally, a comparison of the SMM 
model to the SPI data points yields a x 2 value of about 
8.6 for 6 degrees of freedom (the model amplitude was not 
fitted). 

4. Summary and conclusion 

We studied the sky distribution of Ps continuum emis- 
sion. The extent and shape of extended emission (con- 
sisting of Ps continuum and Galactic continuum) in the 
Ps analysis intervals between 410 and 5 00 keV agree well 
with t hose found for the 511 keV line bv lKnodlseder et alJ 
( 2005). Our bulge parameters are similar to, albeit slig htly 
smaller than, those obtained by IKinzer et alJ (|200 lh for 
the positron annihilation radiation; our Ps continuum map 
is b ulge domina t ed as are the OSSE maps obtained by 
e.g. IMilne et all l|2001bl |2002|) . The emission appears to 
be diffuse, we do not yet find evidence for significant 
emission in the Ps analysis intervals from point sources 
in the GC region. Modelling the extended emission with 
a Gaussian bulge of 8° FWHM located at the GC, and 
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Fig. 4. Comparison of the SPI result for the diffuse emis- 
sion from our Galaxy obtained with a spatial model con- 
sisting of an 8° FWHM Gaussian bulge and a CO disk 
(data points) with the best fit model for the e mission from 
the in ner Galaxy obtained with SMM data bv lHarris et alJ 
l)l990|) (the solid and dashed lines represent their best fit 
model and its total uncertainty). Details are given in the 
text. 



a disk component tracing the Galactic CO distribution, 
we derive spectra and fluxes that are consistent with 



2005; Knodlseder et al. 2005), OSSE 



(Purccll ct al. 


1996; 







observations. 

Emission between 410-500 keV from around the GC 
region is dominated by the Ps continuum. Hence the sky 
distributi on of Ps continuum emi ssion and 511 keV line 
emission ( Knodl seder et alJ l2005|) are consistent, which 
is not surprising considering that most positrons annihi- 
late through formation of Ps. This agreement strength- 
ens our conclusions regarding the origin of Galactic 
positrons based on ob servations of the 511 keV line 
i|Knodlseder et all [2005) . In particular, our results sug- 
gest that an old stellar population, such as Type la su- 
pernovae, classical novae, or low-mass X-ray binaries, is 
the main source of Galactic positrons. Light dark matter 
annihilation and/or decay is another intriguing but spec- 
ulative positron source. 
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